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■ Abstract 

In a relativistic Brueckner-Hartree-Fock calculation we include explicit 
^| negative-energy states in the two-body propagator. This is achieved by using 

\ the Gross spectator-equation, modified by medium effects. Qualitatively our 

qq ' results compare well with other RBHF calculations. In some details signifi- 

cant differences occur, e.g, our equation of state is stiffer and the momentum 
dependence of the self-energy components is stronger than found in a reference 
CN ! calculation without intermediate negative energy states. 
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I. INTRODUCTION 



The relativistic Brueckner-Hartree-Fock model is very succesful in describing equi- 
librium properties of nuclear matter. Where the conventional Brueckner-Hartree-Fock model 
is not able to reproduce the empirical saturation point without the introduction of explicit 
three-body forces, the relativistic extension finds a saturation point very close to the ex- 
perimental value. Going beyond the relativistic Brueckner scheme and including hole-hole 
propagation to all orders, one finds only a small shift in the saturation point towards the 
empirical one ||. Including hole-hole propagation is essential for obtaining a 'conserving' 
model, i.e. a model that preserves conservation laws and is thermodynamically consistent 
@ . In this extended relativistic Brueckner model other observables also compare favourably 
with experiment: we found a compressibility in the range 250-300 MeV, a depletion of the 
Fermi-sea of 11 — 13% and a quasi-particle strength at the Fermi-surface in the order of 0.7. 

The essential ingredient in the RBHF model is the decomposition of the self-energy in 
a scalar and vector part. While the mean-field, in essence the sum of both parts, is in the 
order of ~ 50 MeV, the components are a few hundred MeV. This leads to a small effective 
mass, which again results in a squeezing of the scalar (attractive) part of the interaction. 
This additional density dependent, purely relativistic, repulsive effect is responsible for the 
shift in the saturation point as compared to the non-relativistic Brueckner model. The use 
of relativistic kinematics has only a minor impact. 

One problem intimately related to relativity that has still to be addressed is the effect 
of negative energy states. These are affected by the medium as well, and the question 
to what extent this changes the results is largely unexplored territory. If one takes the 
effective Lagrangian approach on which the RBHF model is based seriously, one should 
attempt to develop some kind of renormalization scheme and chart the effects of the medium 
on this. This has been done in the relatively simple framework of the Walecka model 
in the loop expansion and was only successful in the sense that counterterms etc. can be 
calculated and used on a mean-field level in Hartree and Hartree-Fock approaches on a purely 
phenomenological level 0. However, this loop expansion shows no signs of converging: the 
coupling constants change with every additonal order included || . In a way this failure can 
be expected because the energy scales involved in any renormalization scheme will be in the 
order of a few GeV. On this scale one starts to be sensitive to the substructure of the nucleon 
and therefore to the way how one suppresses this by e.g. a formfactor. Consequently, one 
should not try to renormalize models that are based on point-like nucleons. 

In this work we take a less fundamental approach and avoid the renormalization prob- 
lem by using an effective theory in the vacuum. We employ a model that includes prop- 
agation of negative energy states in a well-defined manner. This is implemented by using 
the so-called spectator equation |J. In this particular three-dimensional reduction of the 
full Bethe-Salpeter equation, negative-energy propagation is explicitly included. In Ref. [T^ 
an interaction is presented that fits the nucleon-nucleon data well using this equation. By 
fitting the data in vacuum one renormalizes the theory 'effectively'. This of course means 
that one loses the dynamics of the renormalization scheme in the medium. However, one 
still has the dynamics of the propagation of the intermediate negative energy state and it is 
interesting to investigate to what extent this affects the conclusions of the standard RBHF 
model. 



2 



We want to stress though that before one can generalize these conclusions more work 
needs to be done. The three-dimensional reduction scheme we use in this work is not unique, 
and therefore any effect will be, at first, specific only to this particular reduction scheme. 
To arrive at more general statements one needs to investigate the effects of negative-energy 
propagation in other reduction schemes. However, before one can do this one needs models 
that fit nucleon-nucleon data including negative energy propagation. We are aware of only 
two efforts in this direction. One is the work of the Utrecht group |TT[ , the other the approach 



of Gross et al. |T(J which we use. 

This paper is organized as follows: we first review in section |TJ the theoretical basis of 
our model. In section [ITT] we will present results for various observables, before we present 
some concluding remarks in the final section. 



II. THEORETICAL FRAMEWORK 

To avoid the complexity of solving the full four- dimensional Bethe-Salpeter equation one 
usually resorts to a three-dimensional reduction scheme. By approximating the propaga- 
tor with a form that fixes the energy variable while preserving the elastic cut, the four- 
dimensional integration is reduced to a three-dimensional one over the spatial momenta. 
Even with the requirement that one retains covariance, there is no unique way to perform 
this reduction. Over the years several specific schemes have been widely used, most notably 
the Blankenbecler-Sugar (BbS) ]I2[ and Thompson ]13[ reductions. In both these equations 



the intermediate nucleons are symmetrically off-shell and no energy is transferred via the 
interaction, i.e., no mesonretardation is present. Another point worth noting is that these 
equations, by construction, only propagate positive-energy intermediate states. Therefore, 
the main difference between these equations and the non-relativistic scattering equation is 
the presence of relativistic kinematics, so-called 'minimal relativity'. RBHF calculations 
based on both equations have been performed by many groups, we mention the pioneering 
work of Anastasio et al. and the calculations by Machleidt 0], Horowitz and Serot |2j and 
the Groningen group ||||. 

A reduction scheme that in some respects goes beyond this is the Gross-reduction, leading 
to the spectator equation |J. In this approach one particle is put on-shell, while in the 
propagation of the off-shell particle both the positive- and negative-energy contributions are 
retained. Also some mesonretardation is present. 

With regard to three-dimensional reductions one should keep in mind that it is not a 
priori clear that the Bethe-Salpeter equation using only the One-Boson-Exchange (OBE) 
interaction as an approximation of the full irreducible kernel is superior to a well chosen 
three-dimensional reduction in describing the full scattering process. For scalar theories it 
is proven that the spectator equation effectively includes the cancellation between direct 
and crossed box diagrams. The spectator equation then gives a better approximation to 
the full sum of ladder and crossed-ladder diagrams than the iterated OBE diagrams of the 
Bethe-Salpeter equation [[|. 

In the spectator approach the 2-body propagator in the cm. frame reads 



,r 2 
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The indices 1,2 identify the particle, in this propagator particle 1 is on-shell. Also, p is the 
cm. 3-momentum, s the total invariant mass and A ± are the usual projection operators 
on positive /negative energy states. It is easily verified that this propagator is covariant. 
This also holds separately for both the positive and negative energy part. Note that we can 
obtain the Thompson equation by ignoring the negative energy part of the propagator and 
the meson retardation. 

In Ref. [TO an OBE model based on the spectator-equation is presented. Using only four 
mesons (ir, p, to and cr) a fit comparable in quality to the Bonn potential |4]] was obtained. 
Several potentials are given, differing in a) taking the pion-nucleon interaction as purely 
pseudo-vector or with an admixture of pseudo-scalar and in b) the treatment of spurious 
poles in the meson propagators due to particular form of meson-retardation in the spectator 
equation. The difference between a rather simple and a more elaborate scheme for the 
treatment of the meson poles were shown to be relatively small. Since our calculations are 
numerically already very involved we opt for the simple scheme and take the IA model of 
Ref. |10j for our calculations in the nuclear medium. 



In isotropic symmetric nuclear matter the self-energy of the nucleon has in the spinor 
representation several components 

E(p) = £ 8 (p) -7°S°(p) +P-1^{ V ). (2.2) 

This decomposition of the Dirac structure allows us formally to define effective (in principle 
momentum and energy dependent) quantities 

m* = m + S s (p) 

p*=p[l + £»]. (2.3) 

Previous calculations have shown that the momentum dependence of the self-energy com- 
ponents is relatively weak. We follow the standard procedure to neglect the momentum de- 
pendence of the self-energies in the intermediate propagators by using the values at p = pf. 
Furthermore, from a technical viewpoint it is advantageous to 'divide out' the E v component. 
One then obtains a renormalized effective mass: 

< = ^±gM , (2.4) 

and corresponding expressions are obtained for momentum and vector self-energies. This 
effective mass will be used in the practical calculation, where confusion might arise we will 
refer to this particular definition as the self-consistent effective mass. With this definition 
we find an effective positive-energy spinor: 



1 



= f^?H I L£ Ur. (2.5) 



2m* 



E* + m* 



Here E* = yp 2 + m* 2 is the on-shell energy of the spinor. As we will see later the value of 
S v is still small, but not negligible as in earlier calculations and we have to incorporate it in 
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the calculation. One thus finds for the dressed one-body propagator |14| (note that in this 
expression p$ is also renormalized by 1 + X v (py)) 



1 m 

g{p) = - 



A*+(p)0(p - p f ) + A* + (p)Q(p f - p) A*~(-p) 



pi- E*+ is Pq — E* — is Pq + E* — is 



(2.6) 



Where A* + is the projector on effective positive-energy states: u*(p)®u*(p), with the effective 
spinors of Eq. Q2.5|) . In line with the Brueckner scheme we neglect the hole (p < pf) part in 
the two-body propagator in the medium. We find 



m* 2 



(1 + h v )E* V vs* - 2E* + is 



JS* — IE 



(2.7) 



In this equation Q PP (P, y/s*, p), Q P (P, vs*, v) are the relativistic angle-average Pauli-blocking 
operators H, respectively blocking out two and one momenta below the Fermi-momentum. 
They depend on the total invariant effective mass vs*, the total momentum P and the 
blocked momentum p. Both momenta in the particle-particle part of the propagator are 
required to have momenta above pf, in the particle- antiparticle part only the momentum of 
the (positive energy) particle is Pauli blocked. 

Any incoming, outgoing or intermediate state can either be particle-particle (denoted by 
'+') or particle-antiparticle (denoted by '-') character. Splitting the two-body propagator in 
two contributions using this nomenclature we find a set of coupled equations for the various 
components of the T-matrix: 

^dir + = Klir + + Klir + 5 , l"^dir + + KxA^^exch + Klir 9 1 ^dir + + Kxch^ ^exch, 

^dii + = Klir + + Klir^l^dir*" + Kx^^exch + Kiir 9l ^dir + + Kxch^ ^exch" (2-8) 



The terminology dir /exch arises from which particle is on-shell (see also Ref. |10|): 



is an interaction where both on-shell particles are on the same-vertex, in each vertex 
has one on-shell particle. In the equation above the outgoing upper particle is always on- 
shell. Which particle is on-shell in the intermediate state is indicated by the index in the 
two-body propagator. The expressions for T^ h have a similar structure. In Ref. [JTO]] it is 



shown that the resulting scattering amplitude is properly anti-symmetrized. By using the 
form of the positive /negative energy projection operators (A* + = u*(p) <S> u*(p)) the above 
set of equations can be expressed in momentum spin-space. In the cm. frame we then 
have, e.g., V++ = u* 3 (p')u{(~p')Vul(p)u* 2 (-p) and T+T = u* 3 (p')v*(-p')Tul(p)u^(-p) and 
so on. The T-matrix is calculated in this representation. We like to stress again that V is 
sandwiched between effective spinors. Due to the low effective nucleon mass this leads to 
a supression of the attractive scalar part and an enhancement of the repulsive vector part 
of the interaction. This addional effective (density dependent) repulsion is a consequence 
of the relativistic treatment and not present in non-relativistic Brueckner models, it is the 
reason for the empirical succes of the relativistic Brueckner models. 

Due to the kinematics introduced by the specific form of the spectator equation, V exch 
can develop a spurious meson-production pole. To counter this we use scheme I of Ref. |TU 



It amounts to replacing the 'exchange' kinematics by 'direct' kinematics in the region where 
the spurious pole can occur. This can be achieved preserving the continuity of V^£. 
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With the T-matrix we calculate the self-energy in the Bruekner approximation: 



dp 



Re 



g<{p). (2.9) 



With T\ the antisymmetrized T-matrix in spinor-space, the '+' indicates that we only take 
the part which is analytic in the upper complex y/s plane. In order to obtain the Brueckner 
contribution to the self-energy we take for the correlation propagator g < (p) only the positive 
energy part 

777* 

g<(p) = 2n— A*+(p)5(p* -E* p ). (2.10) 

By using this expression for the self-energy we eliminate contributions from the Dirac-sea 
and avoid the problems due to the renormalization of these. 

As is clear from Eq. (|2.7|) the full T-matrix has poles in the upper and lower complex 
y/s* plane. The poles in the upper halfplane are generated by the particle-antiparticle 
propagating part, the poles in the lower halfplane by the particle-particle propagating part. 
These poles are well separated as a function of y/s* and we can split the imaginary part of 
T in parts with poles only in the upper or lower half plane 

Im(T) = Im(T + ) + Im(T~), (2.11) 

where T +,_ are analytic in the upper and lower half-plane respectively. This allows us to 
apply Cauchy's theorem for both contributions separately by either closing the contour over 
the upper or lower halfplane. From Eq. (|2.7|) it follows that Im(T~) vanishes on the real 
axis, apart from the point y/s* = 0. This is a highly singular point, we avoid it by using the 
contour depicted in Fig. [1]. Straightforward application of Cauchy's theorem gives for real 



Re 









T-( 


\/s*~) 


=--L 






7T JC 



Im ^£ 1 . (2.12) 



The integrand is zero everywhere except on the halfarc of the contour. There we have to 

calculate the T-matrix with complex values of y/s* in the propagator. We now calculate 
Re(T+) by 

Re(T + ) = Re(T) - Re(T~). (2. 13) 



We tested this method by setting the particle-particle propagating part of the 2-body prop- 
agator of Eq. (|2.7|) equal to zero. Then Re(T + ) should be equal to the Born contribution. 
As it turns out, we need only 6 integration points on the arc to obtain very good accuracy. 
Another test is setting the particle-antiparticle propagating part equal to zero, than Re(T~) 
calculated via Eq. ( |2.12| ) should vanish as well, which is again very well fulfilled. 

In Eq. (|2.9|) the T-matrix is needed in spinor representation. However, in the practical 
calculation we have the T-matrix in spin representation. To extract the spinor representation 
of the self-energy several (approximate) methods have been used in the past. Assuming 
momentum-independent self-energies one can extract the self-energy components by the 
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momentum dependence of the single particle energy j|. Another possibility is to project 
the T-matrix elements in the positive-energy sector of the spin space on Lorentz invariant 
amplitudes as is done e.g. in Refs. Both methods are not exact, for the projection 

method to be unique one needs to include also the negative-energy sector of the spin space 



n|. In the present calculation we want to decompose the self-energy unambiguously. This 
is achieved by calculating the full set of quadratic forms u*(p)Ti(p)u*(p), M*(p)S(p)t>*(p) = 
v* (p)J2(p)u* (p) and v* (p)Y,(p)v* (p) . Working out the straightforward Dirac algebra in a 



similar fashion as in Refs. |1,15| we find 



S s (p) = \ iu*(pMp)u*(p)+v*(p)Z(p)v*(p)] 

= -^K [u*m{p)u*{p) - v*m(p)v*m + -^r g^k^, 

2m* m* 

S v (p) = --L {u*(p)Z(p)u*(p) - v*(p)Z(p)v*(p)] + 3Lv*(p)i:(p)u*(p). (2.14) 
2m* m*p 



Finally, after working out the integration in Eq. ( |2.9| ) and transforming to cm. defining 
variables (s*,P) we have to calculate (note that T + (s*,P) transforms covariantly) 

"^miri ^min ^ P p q 

(2.15) 

In this equation w r and y r can be either a positive or negative-energy spinor, r is the spin- 
index, q is the momentum of the integrated particle in the nuclear matter rest-frame, found 
by applying the inverse transform to the set (s*,P). The integration limits are given by 
^max = P+Pf, Pmm = max(0,p - p f ), s* min = (E* + E* m J 2 - P 2 with q m in = \P - p\ and 
s max = 4-E* 2 — P 2 . The trace in the integrand is over the spin- index s of the integrated 
particle. 

Before we turn to the results, one remaining point of the model has to be addressed. In 
the OBE interaction of Ref. [K| on top of the common mesonic formfactors also a nucleonic 
formfactor, depending on the off-shellness of the nucleon, is introduced. It suppresses large 
off-shell energies and momenta which probe sub-nucleonic structures and reads 

, / 9n 2(A^ - ml) 

k (P) = (K 2 2?, ( K2 2V 2 - 16 

(A2 - p 2 ) + (A2 - ml) 

The question is how to treat this in the medium. Since the formfactor lacks a microscopic 
foundation (apart from the fact that it can be related to an energy-dependent self-energy), 
we have little guidance on how it changes in the medium. Since the on-shell relation in the 
medium reads p* 2 = m* 2 it seems appropriate to use 



2(Af - m. 

(A* 2 - p* 2 ) + (A* n 2 - ml, 
where we scale A* with the effective mass. This idea is slightly generalized by using 



A;= a( <- m " )+m "A n . (2.18) 
m n 
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Taking a equal to zero leaves A* equal to the vacuum value, a = 1 fully rescales A* with 
the nucleon effective mass. 

This completes the theoretical framework of our model. Given a density we now can 
iterate around the self-energy evaluated at the Fermi-surface until a self-consistent solution 
is found. With these self-consistent values we can calculate relevant observables. 



III. RESULTS 

In this section we present results for our model with explicit intermediate negative energy 
states. We compare our results with those obtained using the Thompson reduced version 
of the Bonn-C potential |3J]. The latter model only contains positive-energy intermediate 
states. However, both models do fit the nucleon-nucleon scattering data in the vacuum 
equally well. As we will see, the differences in the underlying dynamical description will 
lead to differences in the results. Of course, the potentials of the models are different as 
well, but this is a consequence of the different dynamics from the requirement that one fits 
the nucleon-nucleon scattering data. 

We first calculate the equation of state, the results are presented in Fig. |[ The full line 
obtained with the cutoff-renormalization parameter a = 1, the dashed line with a = 0.95. 
We see that the dependence on this parameter is (unfortunately) rather large. Moreover, 
lowering a, and thereby letting A n approach its vacuum value, makes the system more 
attractive. If the system gets too attractive, which happens at normal nuclear density 
for a ~ 0.8, it develops an instability which precludes a calculation of the T-matrix. As 
a reference we also plot the result obtained with the Thompson-reduction version of the 
Bonn-C potential ||. In this calculation we used the uniquely determined self-energies and 
did include the (1 + S v ) factor in the 2-body propagator cf. Eq. (|2.7|). One observes that 
the result from the spectator equation is much stiffer than the Bonn-C result. This is partly 
expressed by the compressibility at saturation, 180 MeV for the Bonn-C result and 250 MeV 
for the a = 0.95 result. However, one should be careful with comparing compressibilities 
calculated for different saturation points. The figure does seem to indicate that for higher 
densities the equation of state of the spectator equation is stiffer. A probable explanation 
for this is due to the low effective mass. The particle-antiparticle propagating part provides 
a large part of the repulsion |IIJ. With a lower effective mass, the incoming \/~s* becomes 
smaller as well and the T-matrix equation scales in a sense: the pole of the antiparticle- 
particle part of the propagator comes closer to the pole of the particle-particle part. This 
leads to a larger influence of the particle-antiparticle propagating part and therefore to more 
repulsion. This is an additional repulsive effect of the low effective nucleon mass, on top of 
the conventional relativistic effect of squeezing the attractive scalar part of the interaction. 

We want to stress that this additional repulsion is a direct consequence of the particular 
three-dimensional reduction we use. As pointed out in Ref. [jTOj in fitting the nucleon- 



nucleon scattering data the dynamical repulsion of the particle-antiparticle intermediate 
states results in a relatively low coupling constant of the (repulsive) tu-meson as compared 
to other potentials. Now, only part of the total repulsion needed to fit the data has to come 
from the w-meson. In the Bonn-C on has g^/47r = 20 while in our OBE we use g^Z^ir = 9.85. 
The 'difference' between these values is now included dynamically in the model, and the 
medium effects on this result in differences as we pointed out above. Note that a lower 
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g u brings it in better agreement with the values obtained in pion-nucleon photoproduction 



models, see e.g. Ref. fll6|JT7 



In Fig. |3| we show the density dependence of the effective mass and self-energy compo- 
nents. Our findings are in line with conventional Dirac-Brueckner results. The value for the 
effective mass is similar to what is found in other calculations, we also find the usual linear 
dependence on the Fermi- momentum. The same is valid for E s and S°. The values for these 
are a bit smaller than is usually found, this is mainly due to the fact that we determine the 
self-energy components uniquely. We also find that the isospin 1 = 1 channel contributes 
the most to the respective components of the self-energy. At normal nuclear matter density 
S s = -237 MeV, E° = -197 MeV in the 1 = 1 channel, while in the / = channel the 
respective values are S s = —49 MeV and E° = — 10 MeV, although the contribution of the 
channels to the mean-field (w*Eii*) at the Fermi-surface is comparable: —25 MeV and —6 for 
the I = 0, 1 channels respectively. In an earlier calculation we found the same when using 
the projection method in extracting the self-energies 0,0. A notable difference with other 
calculations is the value of E v : we find a positive and rather large value. The positive value 
is a consequence of the proper, unique, decomposition of the self-energy, similar results are 
found in Refs. fl5Ul4| and also in our recalculation of the Bonn-C results. The absolute size 
of E v is a particular feature of our calculation. It makes the inclusion of the 1 + E v factor 
in the 2-body propagator Eq. (|2~Tj ) necessary. 

It is worthwhile to note that the use of a non-unique method of extracting the self-energy 
components like e.g. the projection method of Refs. ||,|3| and the momentum- dependence 
method of Ref. || leads to differences which, up to order C(p 4 /m* 4 ), only depend on one 
parameter. Two different decompositions of the self-energy E, E' need to give the same 
value for the mean-field: u* (p)T,(p)u* (p) = u*{p)Y/(p)u*(p). From this is follows that up to 
0{j) 4 /m* A ) Yf 3 = E s — a, E' = E — a and T! v = T, v — a/2m*. The 'ambiguity' coefficient 
a depends of course on density and momentum. For example, the momentum- dependence 
method sets E^ = 0, which implies that a = 2m*E„, with E„ the exact component. The 
projection method tends to give negative values for H' v , whereas the exact value is positive. 
This gives an even larger value for the 'error' coefficient a. Our unique decomposition now 
allows us to calculate the mean-field effective coupling constants in a more proper way than 



we did in Ref. ||18|| . In Fig. f| we show the effective isoscalar scalar and vector coupling 
constants. These are defined by 

^X 2 = 1 E' n jp f ) + Elip f ) fg^\ 2 = 1 S°(p,) + YP p (p f ) 
mj ' 2 p s ' VmJ ^ 2 p v ' 1 ' ' 

where p s and p v are the usual scalar and vector densities from relativistic mean-field theory 
T9fl . Comparing the Bonn-C result with those obtained in Ref. |18| using the projection 



method we find that the exact values are smaller and show less dependence on the density. 
The effective scalar coupling only decreases slightly as a function of density, the effective 
vector coupling even increases slightly with density. Both effective couplings are still linear 
in the Fermi-momentum. The appreciable differences between the projection method and 
the exact results show that for detailed results one can not rely on an approximative method 
and has to calculate the unique decomposition. The result of the present calculation with 
the spectator equation shows a much stronger density dependence, also linear in the Fermi- 
momentum. The most striking feature is that both couplings increase with increasing density 
This again reflects the additional dynamics due to negative energy intermediate states. 
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The momentum dependence of the self-energy components is shown in Fig. [| for normal 
nuclear matter density, corresponding to pf = 0.27 GeV/c. In this calculation we take 
a = 1, note that we now present pH v (p). As a reference we give in the upper panel the 
results for the Bonn-C potential and in the lower panel the results of the present work with 
the spectator equation. For the Bonn-C results the conventional assumption of the weak 
momentum dependence of the self-energy, which we used as well, seems to be justified, at 
least for momentum within and slightly above the Fermi-sea. For example, S s increases from 
—320 MeV for momenta deep in the Fermi-sea to —297 MeV for p = pj and for the highest 
momentum give we find —234 MeV. The results obtained with the spectator equation show 
a more pronounced momentum dependence, up to a point where one might start to question 
the usual assumption of weak momentum-dependence of the self-energy components. We 
get for small momenta X s = —344 MeV, at the Fermi surface we have —286 MeV, while for 
the largest momentum we obtain —186 MeV. Again we see that, compared with the Bonn-C 
results, the contribution from E v is significantly larger for the present results obtained with 
the spectator equation. 

IV. CONCLUSIONS 

We presented a relativistic Brueckner calculation using a specific three-dimensional re- 
duction of the Bethe-Salpeter equation that includes explicit negative energy-states in the 
intermediate 2-body propagator. We employ a potential that fits the nucleon-nucleon data 
well using this equation as presented in Ref. [fTOf] . The additional dynamics of the included 
negative-energy states lead to an effective repulsion, expressed by a relatively low tu-meson 
coupling constant. To obtain the Brueckner contribution to the self-energy we need to 
separate out the part of the T-matrix which is analytic in the upper half-plane, i.e. the 
part generated by propagation of particle-particle states. In our calculation we uniquely 
decompose the self-energy in its Dirac structure, removing the ambiguity in determining 
the components present in other calculations We find a binding energy curve that 

compares well qualitatively with previous calculations. This curve does depend significantly 
on how we treat the off-shell form-factor of the nucleon. There are, however, differences 
in important details. Our equation of state if stiffer than found in other Dirac-Brueckner 
calculations. We argued that this is probably an additional effect of the low effective mass. 
This results in the particle-antiparticle propagation, which is repulsive, to become more im- 
portant. The extended dynamics introduced by the particle-antiparticle propagation, which 
is a direct consequence of the particular 3-dimensional reduction of the full Bethe-Salpeter 
equation we use, also expresses itself in the medium. Also, and probably related to this, we 
find the momentum dependence of the self-energy to be relatively large. It is nearly twice 
of what we find when using the Thompson-reduction version of the Bonn-C potential. 

One troubling feature of the model is the medium effect of the off-shell nucleon formfac- 
tor. There is no microscopic description of this formfactor, and hence we do not know its 
medium dependence. We used a simple renormalization of the cutoff and found reasonable 
results. However, the results do depend on the way how we treat the cutoff in the medium. 
This reflects the fact that in our effective theory we are sensitive on how we suppress ener- 
gies that probe sub-nuclear structures. We can only fix this at one density, taken to be the 
vacuum by reproducing the nucleon-nucleon scattering data. However, one then loses the 



10 



dynamical description of the procedure in the medium. This question needs to be addressed 
in future studies. 
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FIGURES 



FIG. 1. The contour used to extract the value of Re(T ) as explained in the text. 

FIG. 2. Binding energies for various values of the rescaling parameter a: the solid line is the 
result with a = 1, the dashed line is calculated with a = 0.95. For reference the result for the 
Bonn-C (Thompson reduction) using the unique decomposition of the self-energy is given by the 
dotted line. 

FIG. 3. Density dependence of the self-consistent effective mass (solid line), T, s (pf) (dashed 
line), T,°(j)f) (dotted line) and £ v (pj) (dash-dotted line). Note the relatively large and positive 
value of £ v (p/). 

FIG. 4. The effective isoscalar scalar and vector coupling constants as defined in the text. 
The solid line is the result obtained with the spectator equation, the dashed line is found using the 
Thompson version of the Bonn-C potential. 



FIG. 5. Momentum dependence of the self-energy calculated at normal nuclear matter density, 
corresponding to pf = 0.27GeV/c. In the upper panel the results for the Bonn-C potential are 
shown, in the lower panel the results for the spectator equation with negative energy intermediate 
states. In both panels the solid line stands for T, s (p), the dashed line for and the dotted line 

for pI7(p). 
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